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ABSTRACT 

We present an analysis of the X-ray spectrum and long-term variability of the nearby dwarf starburst galaxy 
Henize 2-10. Recent observations suggest that this galaxy hosts an actively accreting black hole with mass 
~10 6 Mg. The presence of an AGN in a low-mass starburst galaxy marks a new environment for active galactic 
nuclei (AGNs), with implications for the processes by which “seed” black holes may form in the early Universe. 

In this paper, we analyze four epochs of X-ray observations of Henize 2-10, to characterize the long-term 
behavior of its hard nuclear emission. We analyze observations with Chandra from 2001 and XMM-Newton 
from 2004 and 2011, as well as an earlier, less sensitive observation with ASCA from 1997. Based on detailed 
analysis of the source and background, we find that the hard (2-10 keV) flux of the putative AGN has decreased 
by approximately an order of magnitude between the 2001 Chandra observation and exposures with XMM- 
Newton in 2004 and 2011. The observed variability confirms that the emission is due to a single source. It 
is unlikely that the variable flux is due to a supernova or ultraluminous X-ray source, based on the observed 
long-term behavior of the X-ray and radio emission, while the observed X-ray variability is consistent with the 
behavior of well-studied AGNs. 

Subject headings: galaxies: active — galaxies: dwarf — galaxies: evolution — galaxies: individual (Henize 
2-10) — X-rays: galaxies 


1. INTRODUCTION 


The nearby (~9 Mpc) dwarf starburst galaxy Henize 2-10, 
exhibits intense star formation (e.g., Allen et al.|1976|l, while 
in the center of the galaxy, an X-ray point source (|Kobulnicky 
& Martin 2010) and relatively luminous radio point source 


( Kobulnicky & Johnson) 1999 |Johnson & Kobu l nicky|2003 i 
were found to be co-spatial, suggesting the existence of an 
accreting low-mass active gala ctic nucleus (AGN) with black 


hole (BH) of mass ~1O 6 M 0 (Reines et al. 


2011 


Reines & 


Deller 2012). This represented the first possible detection of 
an AGN in a dwarf starburst galaxy. Even if a large fraction 
of dwarf galaxies host massive BHs, they are challenging to 
detect as AGNs because the AGN emission is faint and its sig¬ 
natures can be swamped by surrounding star formation (e.g., 
Reines et al. 2013|); X-ray observations can be one of the most 


effective methods for identifying AGN in dwarf galaxies (e.g., 


|Reines et al.||2014~t [Lemons et al.| 2015[ |Secrest et al.|2015) . 

If the existence of an AGN in He 2-10 is confirmed, it would 
serve as one of the best possible analogs for BH and g alaxy 
grow th in the early history of the Universe (e.g., Reines et al. 
[201 l) l. Most bulge-dominated galaxies contain supermassive 
BHs, however, the process by which the orignal “seed” BHs 
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formed remains poorly constrained (e.g., Johnson & Bromm 
[20071 |Volonteri|2010) . 

Currently, the available observational evidence for the cen¬ 
tral compact sources in He 2-10 favors its interpretation as 
a supermassive BH. The majority of its radio emission orig¬ 


inates from a region < 3 pc xl pc in size (Reines & Deller 


2012|l and is consistent with being spatially coincident with 
Ci 


the Chandra hard X-ray point source at the dynamical cen¬ 
ter of the galaxy (Reines et al.] |2011) . Assuming that the ra¬ 
dio and X-ray emission are produced by a BH, a c omparison 
with the BH fundamental plan e (Merlon i et ak |2003 i suggests 
that the mass is ~10 6 M 0 (Reines et al.|2011 1 . Alternatively, 
the X-ray emission could in principle come from an ultralu¬ 
minous X-ray source that is powered by a stellar-mass BH 
( |Roberts|2007| l. Howeve r this cannot account for the oberved 
compact radio flux (e.g. [Middleton ^t al.||2013[ |Wolter et al.| 


2014|l, although we note that previous radio and X-ray ob- 


servations are not strictly simultaneous. We can most likely 
rule out supernova (SN) remnants as the cause of the X-ray 
emission; there are no massive star-forming clusters coinci¬ 
dent with the compact radio emission, rendering this scenario 
somewhat implausible, but not impossible. To more robustly 
constrain the nature of the compact central source in He 2-10 
and to better constrain its mass, it is important to understand 
how its X-ray luminosity varies with time. This is the goal of 
the present paper. 

The original evidence for the AGN in He 2-10 came in part 
from analysis of the spectrally hard, resolved point source in 
the 2001 Chandra observations. Here, we analyze data taken 
from Chandra (2001), XMM-Newton (2004 and 2011), and 
ASCA (1997) to obtain spectra at each epoch and a result¬ 
ing measure of the long-term variability of the hard nuclear 
source. The temporal baseline of the observations is suf¬ 
ficient to probe variability on timescales reasonable for an 
intermediate-mass BH or low-mass AGN, as shown by the 
small known sampling of these rare objects (Jpewangan et al. 
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TABLE 1 

Observation Details 


Instrument 

Detector 

Obs. ID 

Date 

Exposure in ks (clean) 

Net counts 

0.5-2 keV (2-8 keV) 

Chandra 

ACIS-S 

2075 

2001-03-23 

20.0(19.7) 

983(174) 

XMM-Newton 

P n 

0202650101 

2004-05-27 

42.0 (29.3) 

3216(234) 

XMM-Newton 

pn 

0672800101 

2011-05-11 

26.9(17.6) 

1863 (92) 

ASCA 

SIS 

65017000 

1997-11-30 

39.8 (22.4) 

197 (52) 


Note. — The details of each of the observations used in this paper. Observed net counts after 
background subtraction are listed for the 0.5-2 and 2-8 keV bands (see §[2]for details of source 
extraction and background analysis). Exposure times listed are for the total exposure and the 
net exposure after cleaning for flares. We focus on our analysis on the more sensitive Chandra 
and XMM observations, but include the earlier ASCA observation as a check on the baseline flux 
level for the source. 


12008). 

2. DATA REDUCTION AND SPECTRAL ANALYSIS 

In this section we describe the spectral extraction process 
for each observation, with a focus on details of the sophis¬ 
ticated background modeling, which was required for the 
XMM-Newton and ASCA observations. All spectral analyses 
are performed using XSPEC vl2.8.0 ( |Amaud|l996 1 . See Ta¬ 
ble [T] for the details of the observations. We will focus pri- 
mamy on the more sensitive observations taken by Chandra 
in 2001 and XMM in 2004 and 2011, but will also include a 
discussion of the 1997 ASCA observation that provides a (less 
sensitive) baseline measurement of the source flux. Through¬ 
out the paper, uncertainties on X-ray measurements (fluxes, 
luminosities, and spectral parameters) represent 90% confi¬ 
dence intervals. 


2.1. Chandra 

The nuclear X-ray point source in He 2-10 was discov¬ 
ered in the 20 ks ACIS-S Chandra observation of He 2-10 
on 23 March 2001 (Reines et al .|2011| ). The pipeline-reduced 
data from this observation (see Figure []} were obtained from 
the HEASARC public archive and were reduced using the 
analysis tools of CIAO 4.5. Time filtering yielded 19.7 ks 
of cleaned exposure. For Chandra we defined two separate 
source extraction regions: a small region of radius 2.25" at 
the nuclear region of the galaxy was used to measure the flux 
of the hard point source, and a larger region of radius 16.8" 
was used to include the soft diffuse X-ray component. In ei¬ 
ther case, the background region was comprised of an annu¬ 
lus of outer radius 58 arcsec, excluding the larger (diffuse) 
source region and another point source about 45 arcsec away 
from the source. We note that the background level above 
2 keV in the Chandra observation was much less prominent 
than in the XMM and ASCA observations. The spectrum, as 
well as the response and ancilliary files, were extracted us¬ 
ing the specextract command. Due to the high signal- 
to-noise ratio and low background, we used CIAO’s built- 
in background subtraction rather than simultaneously fitting 
source and background spectra, as for the XMM and ASCA 
observations. 

For all the observations presented in this paper, the 
source spectrum is described with a model consisting of the 
power-law component to model t he hard nuclear source and 
optically-thin thermal (VMEKAL;|Mewe et al.|1985[|Liedahl 
et al.|1995| t) component to model the diffuse emission, with 
abundance values fixed to those obtained in |Kobulnicky & 
Martin ( 2010) ) (0.78 for light elements and 0.29 for heavy el¬ 
ements), and allowing the normalization and temperature to 


float. We use a VMEKAL to match the spectral analysis of 
Reines et al. (2011), but note that fitting with an APEC model 


(Smit h et al.|2001| ) has no sign ificant effect on the results fo r 
the hard component. Following Kobulnicky & Martin (2010), 
we include Galactic absorption An,Gai = 5 x 10 zu cm - ^ on 
all components, and local absorption for the VMEKAL com¬ 
ponent (Ah. Diffuse = 9.7 x 10 20 cm -2 ). The absorption on 
the power law component (Ah,N uclear) was allowed to float. 
Absorption is computed using the tbabs model (Wilms et al. 


2000|. In XSPEC notation, the source model is given by: 


Source = tbabsGai(tbabsDiffuse*VMEKAL + tbabsN U ciear*powerlaw) 

( 1 ) 

We first fit this model to the spectrum from the nuclear 
(2.25" radius) source region, to obtain the strongest possible 
constraint on the emission from the unresolved hard compo¬ 
nent. With relatively few counts at high energies we obtain 
poor constraints on the hard X-ray photon index, so this is 
fixed to the canonical AGN value of T = 1.8. (Repeating 
the fit for values of intrinsic T varying within a range typi¬ 
cal for AGN, 1.4 < T < 2.2 (]Tozzi et al.] [2006) , produces 
no significant change in the unabsorbed flux.) The fit yields 
Ah = (4.61±^) x 10 22 cm -2 , indicating substantial ab¬ 
sorption. We next fit the same model to the spectrum from the 
extended source region (a circle centered on the nuclear re¬ 
gion with 16.8" radius), but fix the Ah value on the power law 
component to that obtained for the nuclear spectrum. We ob¬ 
tain a consistent and nearly identical flux for the hard power- 
law component between the extended and nuclear regions, 
but a substantially brighter diffuse (VMEKAL) component in 
the extended region, with best-fit kx = 0.65 ± 0.03 keV. 

This confirms that for an even larger extraction region such 
as those used for XMM and ASCA, the hard spectral compo¬ 
nent can be associated with the compact nuclear source. The 
best-fit fluxes and spectral parameters are given in Table [2] 

We quote intrinsic (unabsorbed) fluxes and luminosities for 
the hard nuclear component, and observed (absorbed) fluxes 
for the soft diffuse component. For direct comparison with 
the XMM analysis, we have also extracted a Chandra spec¬ 
trum with somewhat larger radius of 36", corresponding to the 
XMM source region described in the next section (Figure IT]). 

Using this larger source region has no significant effect on the 
spectral parameters. 


2.2. XMM-Newton 

We use two subsequent observations of He 2-10 by XMM 
to constrain the long-term variability of the source. The XMM 
observations on 27 May 2004 and 11 May 2011 have expo- 
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Fig. 1.— Chandra images in the soft (0.5-2 keV) and hard (3-8 keV) X-rays; energy ranges are chosen to clearly separate the soft diffuse component from 
the hard compact nuclear emission. The XMM source extraction region, as well as the two Chandra regions are superimposed. The excellent angular resolution 
of Chandra allows for clear imaging of the X-ray morphology. Widespread diffuse emission from star formation is seen in soft X-rays, while the central nuclear 
source is clearly seen at hard X-rays. XMM-Newton has significantly poorer angular resolution, thus a larger extraction region was required in order to include 
sufficient source flux. 


sure times of 42 and 27 ks, respectively. For both XMM ob¬ 
servations, we reduced, cleaned, and extracted spectra from 
all three CCD cameras: pn, MOS1, and MOS2. After spectral 
extraction and analysis, the MOS1 and MOS2 data yielded 
significantly lower signal-to-noise ratios at energies >2 keV 
compared to the pn detector, so that no useful constraints were 
obtained on the hard emission from the nuclear point source. 
In what follows we therefore focus on results from the pn. 

The source extraction regions were 36" in radius, chosen to 
provide a balance between extracting as many counts as pos¬ 
sible from the source and minimizing the background. (As a 
check, we have repeated the analysis using a smaller source 
region of 25" radius, and obtain essentially identical results 
with marginally larger uncertainties.) The Chandra images 
(Figure[l]i show that the extent of both the nuclear and diffuse 
components are substantially smaller than the 36" extraction 
region, such that they can both be considered as approximate 
point sources for the XMM analysis. In both data sets the 
source region was on-axis and did not lie on any chip gaps. 
We extracted pn spectra for counts in the energy range 0.2-15 
keV and with event patterns 0-4. Response files were pro¬ 
duced from the XMM-Newton Current Calibration Files cor¬ 
responding to the time of this observation. The source ARF is 
calculated including a correction for photons falling outside 
the extraction region. This energy-encircled fraction (EEF) 
varies with energy but is ~85% at 5 keV. The source spectrum 
is described with the same model as for Chandra, consisting 
of VMEKAL and power law components (Equation [T]). 

To maximize the number of counts in the background spec¬ 
tra and thus achieve the highest possible signal-to-noise ratio, 
we extracted the background spectrum from a large annulus 
of outer radius 3 arcmin around the source. Based on a num¬ 
ber of trials, the 3 arcmin annulus was determined to provide 
the optimum number of background counts without needing 
to account for the variation in background flux at larger off- 
axis angles. Using the SAS command edetect_chain, 5 
sources in the pn field of view were detected and subsequently 
excluded from the background region. With these sources ex¬ 
cluded, the background region is «23 arcmin 2 in area, or 20 
times larger than the source region. Because the background 
emission is extended in nature, the ARF for the background 
region did not include an EEF correction. 

The pn background spectrum was fitted with two compo¬ 
nents. The instrumental background is modeled by a power 


law continuum, plus Gaussian emission lines caused by fluo¬ 
rescence (A1 K-o at 1.5 keV and CuNi K-a at 8.5 keV ( |Carter| 
& Read|2007]l. Because this instrumental background is pro¬ 
duced internally to the detector and is not affected by the 
mirror response, in modeling the observed counts it was was 
convolved with an RMF but not multiplied by an ARF. The 
background line energies were determined from fitting each 
line individually, and then fixed for the full spectral analysis, 
while the intrinsic line widths are fixed to be consistent with 
zero. 

The sky background component is dominated by the diffuse 
soft cosmic X-ray background (CXB), which can be modeled 
as thin-thermal emission ( |Hickox & Markevitch|2006]|2007j ). 
We used an APEC model for t his c omponent, to match the 
spectral shape obtained by Hickox & Markevitch| (2006) in 
fitting the unresolved CXB spectrurrnn the Chandra Deep 
Fields, fixing kT = 0.17 keV and allowing the normalization 
to float. The hard (> 2 keV) CXB can be described as a power 
law owing to the summed emission from a large number of 
AGN (e.g., Hickox & Markevitch|2006| . We did not include 
this as a separate component here, as the average expected 
numbers of counts is <2% of the instrumental background, 
so its small contribution to the total background flux can be 
effectively accounted for in our modeling of the instrumental 
background. 

The full model for the total observed emission in the XMM 
source region is: 


Data = Source + Instrumental BG + Sky BG. (2) 

The source spectrum is modeled by an absorbed VMEKAL 
and power law (Equation ITb, while the background compo¬ 
nents are modeled, in XSPEC notation, as: 


Instrumental BG = powerlaw + gauss + gauss + gauss 

(3) 

and 

Sky BG = APEC, (4) 

where Instrumental BG is convolved with the RMF only, 
while Source and Sky BG are convolved by the RMF and 
multiplied by the ARF. 

As discussed below, the fluxes of the nuclear power law 
component in the XMM observations are significantly smaller 
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Fig. 2.— X-ray spectra including model fits and residuals, for four X-ray observations of He 2-10. Each component of the model fitted to the data are shown 
as dotted lines: the diffuse VMEKAL component (red) and nuclear power law (blue). For the Chandra observation (a) the background was subtracted before 
spectral fitting, while for the XMM (b), (c) and ASCA (d) observations, the background spectra are fitted by a model(shown by the gray lines) simultaneously 
with fitting of the observed source spectrum. There is a strong and significant detection of the hard nuclear power law component in the 2001 Chandra (a) 
observation (clearly visible as a point source in Figure^. The hard component is significantly weaker in the XMM observations (b), (c) indicating variability by 
approximately an order of magnitude. The hard component is also detected, at lower significance, in the 1997 ASCA spectrum (d). 


than that observed for Chandra. To most accurately extract 
the weak hard X-ray signal from the significant background, 
we modeled the data by simultaneously fitting the source and 
background spectra, using the model given in Equation[2] We 
account for the differences in area between the source and 
background regions by setting the AREAS CAL parameter on 
the background spectrum. The spectral parameters of the in¬ 
strumental and sky backgrounds were fixed to be equal for 
the source and background spectra. The source component is 
also included in the background spectrum, but multiplied by 
a factor 5 x 10 -3 to approximately model the flux scattered 
outside the source region into the background region. (The 
scaling factor accounts for both the energy encircled fraction 
and the relative area of the source and background regions; 
the ultimate fit parameters are insensitive to the precise value 
of this factor.) 

Because the soft (VMEKAL) emission from the source is 
diffuse (with a diameter of «5"or 200 pc), and thus a light¬ 


crossing time much longer than the separation in time be¬ 
tween observations, it should not be observed to vary in our 
data. Therefore, to maximize the statistical power of our 
modeling, we fitted the 2004 and 2011 XMM spectra simul¬ 
taneously, tying the temperatures and normalizations of the 
VMEKAL and APEC components between the two data sets. 
We allowed the the nuclear power-law normalization to float, 
along with the instrumental background parameters (the par¬ 
ticle background in the detectors should not be perfectly con¬ 
stant for the duration of the mission). We thus performed a 
simultaneous fit to four spectra: source and background spec¬ 
tra from each of the 2004 and 2011. 

The results of the spectral fitting are shown in Figurej2]and 
listed in Table [2] We observe a clear decrease in the flux of 
the hard nuclear component between the 2001 Chandra and 
2004 XMM observations. This is demonstrated in Figure RJ 
in which we show the 2004 XMM pn spectrum fitted with the 
2001 Chandra best-fit model, with no model for the XMM 
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Fig. 3.— Spectrum of the 2004 XMM observation, fitted with source spec¬ 
trum from the 2001 Chandra observation. This fit includes no subtraction or 
modeling of the XMM background. This large excess of the model over the 
data at high energies clearly illustrates the dramatic decrease in the flux of the 
hard nuclear component between 2001 and 2004, independent of the methods 
used to account for the XMM background. 


background included. This shows that the hard source flux 
has dropped dramatically from the Chandra level, even before 
accounting for the XMM background. We observe a further, 
less significant decrease between the 2004 and 2011 XMM ob¬ 
servations, while the flux of the diffuse (VMEKAL) compo¬ 
nent is consistent with no variation from the Chandra obser¬ 
vation. We also observe evidence for a decrease in absorption 
on the nuclear component, with Ah,N uclear consistent with 
zero. The best-fit kT of the VMEKAL component for XMM 
is close but not fully statistically consistent with the Chandra 
data (0.58 ^q'q 3 keV compared to 0.65 ± 0.03 keV). Fixing 
this temperature to the Chandra value has a negligible effect 
on the flux of the diffuse component, but decreases the flux of 
the hard nuclear component by ssl5%. This results in an even 
larger observed drop in flux compared to the Chandra data; 
in the following discussion we will conservatively consider 
the smaller change in flux obtained when the VMEKAL kT 
allowed to float for XMM. The flux of the APEC component, 
representing the soft diffuse CXB, corresponds to a 0.5-2 keV 
surface brightness of (2.9±0.5) x 10~ 12 erg cm -2 s -1 deg~ 2 , 
similar to the soft back ground intensity obtained in the Chan¬ 
dra Deep Fields ( Hickox & Markevitchfl2006j >. The overall 
implications of the spectral fitting results are discussed in §[3] 
One potential uncertainty in our spectral analysis arises 
from the implicit assumption that the instrumental and sky 
background in the source region has the same surface bright¬ 
ness to the emission in the background region. It is possible 
that fluctuations in the background could cause this assump¬ 
tion to be invalid, leading to over- or under-subtraction of the 
background, particularly at energies >2 keV where the back¬ 
ground dominates the signal. To directly check the level of 
possible fluctations, we extracted the 2-10 keV counts in 80 
circular regions of 36" radius (equal to the source region) sur¬ 
rounding the source region. We avoid chip gaps and obvious 
bright sources, noting that there is no bright point source de¬ 
tected within 36" in the 2-8 keV Chandra image (Figure [III. 
The 2-10 keV counts in these apertures are approximately 
normally distributed, with mean (dispersion) of 239 (19) and 
138 (15) counts for 2004 and 2011, respectively. We conclude 
that our modeling of the background in the source region is 
dominated by statistical error rather than any systematic un¬ 


certainty due to background fluctuations. 

2.3. ASCA 

As a measurement of the baseline level of flux prior to the 
Chandra observations, we also utilize observations of He 2- 
10 from ASCA in 1997. Due to ASCA’s poor angular resolu¬ 
tion and sensitivity relative to Chandra and XMM , it provides 
relatively weak constraints on the source flux, particularly in 
the hard nuclear component. We will therefore focus our con¬ 
clusions primarily on the Chandra and XMM data, but will 
utilize the ASCA data here as a useful check on our conclu¬ 
sions. 

He 2-10 was observed by ASCA for a total of 39.8 ks on 30 
November 1997. We used Xselect v2.4b to remove time inter¬ 
vals of high background for a net exposure of 22.4 ks, and to 
extract source and background spectra from the cleaned event 
files. We extracted a spectrum in the energy range between 
0.5-10.5 keV, using a source extraction region of 1.45 arcmin 
in radius. This source region was chosen to contain as much 
of the source flux as possible (the energy encircled fraction 
is 50%) while minimizing the contribution from background. 
We extracted a background spectrum from a rectangular annu¬ 
lus of area sa5 times that of the source region, located around 
the source position but excluding the source region. The RMF 
and ARF response files were generated for each chip using 
the commands sisrmg and ascaarf, respectively. Fur¬ 
thermore, the spectra from the two chips SIS0 and SIS 1 were 
added using the HEASOFT command mathpha, while the 
response files were added with area dependent weights u sing 
the commands addarf and addrmf. As discussed in § |2.2[ 
the source ARF includes a correction for the energy encircled 
fraction, while the background ARF does not. 

As with XMM, we modeled the ASCA data by fitting the 
source and background spectra simultaneously, modeling the 
scattered flux in the background region assuming an energy- 
encircled fraction of 0.5. The instrumental background was 
modeled with a power law plus five Gaussian emission lines; 
three were included in the model as narrow fluorescence lines 
stemming from the device itself at 6.5 keV, 7.5 keV, and 8.2 
keV, which are the Fe and Ni K-a lines and the Ni K-/3 linej^] 
Another Gaussian at 3.3 keV of unknown origin was intro¬ 
duced to fit a feature of the background spectrum, while the 
final Gaussian was a broad line at 11 keV that modeled the 
internal background above 7 keV to account for the steepness 
of the power law. 

We utilized the same source and sky background models as 
for XMM (Equations [T] and [4|. Because they cannot be well- 
constrained owing to the poor photon statistics, the surface 
brightness of the APEC component and the VMEKAL tem¬ 
perature are fixed to the values obtained with XMM, and the 
photon index of the power law is again fixed to T = 1.8. 
(Allowing the VMEKAL temperature to float yields a sig¬ 
nificantly higher kT sa 0.9, inconsistent with the XMM and 
Chandra results, but has a negligible effect on the total fluxes 
of the diffuse and nuclear components.) The results of the 
ASCA spectral fitting are shown in Figure [2] and listed in Ta¬ 
ble [2] We obtain a significant detection or both the diffuse 
and nuclear spectral components, although with significantly 
larger uncertainties than in the Chandra and XMM data. In 
contrast to the XMM observations, the fit yields significant nu¬ 
clear absorption consistent with the Chandra value, although 
the precise value of Ah is poorly constrained. 

9 https://heasarc.gsfc.nasa.gov/docs/asca/newsletters/sis_back2.html 
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TABLE 2 

Spectral fitting results 



VMEKAL (Diffuse) 



Power Law (Nuclear) 



Observation 

kT 

H). 5-3 keV 

Ah, Nuclear 

r 

Lb-10 keV 

L 2 -IO keV d 

xl (d.o.f.) 


(keV) 

( 10 —13 erg cm -2 s _1 ) 

( 10 22 cm -2 ) 


( 10 —13 erg cm -2 s -1 ) 

( 10 39 erg s _1 ) 


Chandra h (2001) 

0.65 ± 0.03 

2.06 ± 0.11 

4 61 +1 ' 67 
—1.26 

[ 1 . 8 ] 

O 90+0.73 
^•^o_0.64 

q 1 O+0.71 
°-- LO — 0.62 

1.61 (46) 

XMM C (2004) 

0 5«+°- 02 
U,OS -0.03 

1 90 +0 ' 08 
i.yu_o 05 

< 0.05 

[ 1 . 8 ] 

o.68i° 0 ;°? 

0 66 +0 ' 08 
U.OO_o 07 

1.20 (661) 

XMM C (2011) 

0 58+ 0 - 02 
U,O5 -0.03 

1 qn+0.08 

i.yu_o.o5 

< 0.05 

[ 1 . 8 ] 

0.44tO;09 

0.43 ±0.09 

1.20 (661) 

ASCA (1997) 

[0.65] 

9 cq+0-89 
—1.02 

u -ou_ o.3o 

[ 1 . 8 ] 

1.79+jS 

1 '70 + 1.61 

1 °-0.84 

0.80 (70) 


NOTE. — Best-fit spectral parameters obtained from modeling of the four X-ray spectra. The (unabsorbed) 2-10 keV fluxes were calculated based on the nuclear 
(power law) component, while the (observed) 0.5-3 keV fluxes correspond to the diffuse (VMEKAL) component. Both components are modified by Galactic absorption 
with column density fixed at Ah, Gal = 9 X 10 20 cm -2 , and the VMEKAL component is absorbed by an additional component with column density fixed at 
Ah, Diffuse = 9.7 X 10 20 cm -2 . Parameters in the table that are fixed in the fits are identified with brackets. Uncertainties represent 90% confidence intervals. 
a Luminosity values were calculated assuming a distance of 9 Mpc to He 2-10. 

k Parameters in the Chandra fits for the nuclear and diffuse components are determined from the fit to the nuclear and extended source regions, respectively, as 
described in m 

c The spectra for the two XMM observations are fitted simultaneously, with parameters for the diffuse (MEKAL) component tied between the two observations, as 
described in §|2.2| 


3. RESULTS AND DISCUSSION 

The long-term X-ray light curve of He 2-10, showing varia¬ 
tions in the VMEKAL and power law components, are shown 
in Figure [4] The diffuse component shows no significant 
variability over the four observations, as expected for emis¬ 
sion from a large-scale diffuse plasma. In contrast, it is im¬ 
mediately clear that there is significant variability in the hard 
power law component. The luminosity of the nuclear source 
decreased significantly between the 2001 Chandra and 2004 
XMM observations, and by approximately an order of mag¬ 
nitude between 2001 and the 2011 XMM observation. (We 
note that our overall conclusions are unchanged if we use ob¬ 
served hard fluxes, which differ by roughly 30% from the val¬ 
ues shown in the Table 2 for the Chandra and ASCA observa¬ 
tion and remain constant for XMM due to the decreased levels 
of obscuration.) 

The variation in the light curve of hard spectral component 
He 2-10 over approximately an order of magnitude in Lx 
confirms that this emission is the result of a single object, 
rather than several separate sources. This allows us to per¬ 
form comparisons to other individual astrophysical sources. 
One class of object that can have similar X-ray and luminosi¬ 
ties and amplitudes of variability is SNe (see Dwarkadas & 
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Fig. 4.— The fourteen year light curve of He 2-10 with 90% confidence er¬ 
rors shown. The hard nuclear flux (blue) declines between the 2001 Chandra 
and 2011 XMM observations by nearly an order of magnitude. In contrast, 
the soft diffuse flux (red), remains approximately constant between the four 
observations. 


Gruszko|2012 for a compilation of published SN X-ray light 


curves). However, a SN interpretation for the nuclear source 
is inconsistent with the radio properties. The radio flux of 
He 2-10 has been measured at 5 GHz with the VLA in 1994 
(0.89 ± 0.18 mJy ; |Kobulnicky & Johnson[|l999l l and 2004 
(0.86 ± 0.02 mJy; Reines et al.|201 lj)7implying no significant 
change in the radio flux withTime. Long Baseline Array mea¬ 
surements in 2011 at 1.4 GHz yield a flux of 0.98 ± 0.21 mJy 
in a compact source on sal pc spatial scales; assuming a typ¬ 
ical radio spectral index, this implies that roughly half of the 
observed 5 GHz flux comes from the central compact source 
( Reines & Deller[2012| >. High Sensitivity Array observations 
in 2005 did not detect t he compact source on extremely small 
(~0.1 pc) spatial scales (|Ulvestad et al.|2007), placing a lower 
limit on its spatial extent. ThisTules out out the presence of 
a single very young SN, however the total radio luminosity of 
the source would imply that any single SN must be at most 
decades old fFenech et al.|2010| . 

With these constraints in mind, we test whether the nuclear 
source in He 2-10 is consistent with a (relatively) evolved 
SN explosion, by comparing its (unabsorbed) soft X-ray and 
5 GHz radio light curves to a sample of 7 SNe that have 
both radio and X-ray measurements in the compilations of 
jDwarkadas & Gruszko| ( |20l2i ) and |Weiler et al.| ( |2002| >. As¬ 
suming that the nuclear source in He 2-10 has a relatively 
constant 5 GHz flux between 1994 and 2004, and conserva¬ 
tively assigning all the 5 GHz flux observed with the VLA 
(«0.9 mJy) to the central compact component, we find that 
the ratio of X-ray (0.5-2 keV) to radio (yF v at 5 GHz) for 
He 2-10 is ~5 x 10 3 . This is more than an order of mag¬ 
nitude larger than the typical X-ray to radio flux ratio for X- 
ray-detected SNe at ages >1 year, and 2.5 times larger than 
the most extreme observed values, in SNe 1980K and 1970G. 
Further, the detection of the nuclear hard X-ray component 
in the ASCA observation indicates that the X-ray light curve 
rises or remains constant, and then declines sharply over a few 
years. This behavior is unusual for SNe; of the eight SN X- 


ray light curves in the compilation of Dwarkadas & Gruszko 
( |2012 > that extend beyond ten years, none shows a similar 
sudden, rapid decline on these time scales. We also note that 
the X-ray to radio flux ratio is consistently at least 2-3 orders 
of magnitude smaller than the typical ratio of X-ray to com¬ 
pact radio flux for ULXs, for which the compact sources tend 


to be weak or undetected in the radio (e.g., Middleton et al. 
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2013 Wolter et al. |2014]l. Given the observed X-ray light 
curve and the ratio of X-ray and radio luminosities, we con¬ 
clude that the observations do not favor a SN or ULX origin 
for the nuclear source. 

In contrast, the significant variability of the hard nuclear X- 
ray source in He 2-10 is consistent with its identification as 
an accreting massive BH, in comparison to the X-ray variabil¬ 
ity of known low-mass AGN. The Sdm spiral galaxy NGC 
4395, at a distance of only 4 Mpc, contains a BH of mass 


3.6 x 1Q 5 Mr.'j , whose hard component, as is shown in King 


et al. (20131, can vary by a factor of 5 on a timescale of just 


one day. The nearby edge-on Seyfert 2 galaxy NGC 4945, 
with a BH mass obtained through observations of its H 2 O 
megamaser of a; 1O 6 M 0 , shows intrinsic variability (mea¬ 
sured at >8 keV by RXTE and Swift/BXl) of at least an order 
of magnitude on timescales of day s to weeks (e.g., [Mueller 
et al. ||2004[ IMarinucci et al.| |2012[ >. The AGN in the n earby 


Seyfert 1 galaxy NGC 4051 (Mrm = 1-73 x 1O 6 M 0 ; 


Den¬ 


ney et al. 2009) has been observed to vary in X-ray luminosity 
by more than an order of magnitude over ~ year timescales 


(Uttley et al. |1999 i. This limited survey confirms that X-ray 
variability over a large dynamic range on timescales of years 
is not uncommon among relatively low-mass AGNs. There¬ 
fore, the decreased observed luminosity after the Chandra ob¬ 
servations could either be due to short timescale fluctuations 
occurring precisely at the time of the observation, as in NGC 
4945, or part of a general trend of long timescale variability, 
similar to objects like NGC 4051. 
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